TUM/T39-97-19 



Gluon Polarization in NucleonsQ 



A. Saalfeld, G. Piller and L. Mankiewicz[] 

Institut fur Theoretische Physik, TU Miinchen, Germany 
February 1, 2008 



Abstract: 

In QCD the gauge-invariant gluon polarization AG in a nucleon can be defined either 
in a non-local way as the integral over the Ioffe-time distribution of polarized gluons, or 
in light-cone gauge as the forward matrix element of the local topological current. We 
have investigated both possibilities within the framework of QCD sum rules. Although the 
topological current is built from local fields, we have found that its matrix element retains 
sensitivity to large longitudinal distances. Because QCD sum rules produce artificial 
oscillations of the Ioffe-time distribution of polarized glue at moderate and large light- 
like distances, the calculation of the matrix element of the topological current results in 
a small value of AG(fi 2 ~ 1 GeV ) ~ 0.6 ± 0.2. In a more consistent approach QCD 
sum rules are used to describe the polarized gluon distribution only at small light-like 
distances. Assuming that significant contributions to AG arise only from longitudinal 
length scales not larger than the nucleon size leads to AG(fi 2 ~ 1 GeV 2 ) ~ 2 ± 1. 
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1 Introduction 



Thanks to celebrated factorization theorems |3J hard scattering of highly virtual probes 
from nucleons can be characterized in QCD by universal, process independent, non- 
perturbative distribution functions which contain all relevant information about the long- 
distance dynamics of the target. At twist-2 these are identified in the framework of the 
QCD-improved parton model with scale-dependent quark and gluon light-cone distribu- 
tions. For spin-1/2 targets they include unpolarized quark and gluon distributions and 
their polarized counterparts related to longitudinal and transverse target polarizations. 
Quantum numbers and chiral properties of the hard probe determine which particular set 
of distribution functions can be accessed in a specific process. Until now most informa- 
tion about the nucleon structure has been obtained from deep-inelastic lepton scattering 
experiments. In particular, recent measurements of polarized deep-inelastic scattering at 
CERN [0] and SLAC || have shown that a relatively small fraction of the nucleon spin 
is carried by quarks. This has started an ongoing debate about remaining contributions 
to the nucleon spin [[| which may result from gluon polarization and orbital angular mo- 
mentum. Here especially the polarized gluon distribution AG(u, /i 2 )Q became of interest 
since it turned out to be measurable in future high-energy experiments, e.g. charm and 
direct photon production || || . 

At present the only available information about the magnitude of the total gluon 
polarization 

AGV) = / duAG{u,fi 2 ) (1) 



results from the analysis of scaling violations in the polarized structure function g±(u, Q 2 ). 
Here a relatively large gluon polarization has been found even at a low normalization scale 
AG(/i 2 = lGeV 2 ) = 1.6 ± 0.9 0. The main objective of the present paper is to discuss a 
framework in which one may determine AG(/i 2 ) using presently known non-perturbative 
methods. 

The link between parton model ideas and QCD is provided by the Operator Product 
Expansion (OPE) |T|. It allows to relate moments of parton distributions to matrix ele- 
ments of local operators with appropriate quantum numbers^. The latter are computable 
either in some approximate way using e.g. model descriptions of hadronic structure M, Rf 



or QCD sum rules |]10| , or at least in principle, through lattice simulations |[L1 |. 

An alternative, but completely equivalent picture views twist-2 parton distribution 
functions as normalized Fourier transforms of nucleon matrix elements of non-local QCD 
operators [JT2J, constructed as gauge- invariant overlap of two quark or gluon fields sep- 
arated by a light-like distance. Since the first few coefficients of the Taylor expansions 
of these non-local matrix elements around the origin are given by matrix elements of 



throughout this paper we denote the Bjorken variable xb 3 by u to avoid confusion with the space-time 
variable x. 

2 We define the n-th moment of a distribution F(u) as T TL [F] = J Q duu n ~ 1 F{u). 
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local, low- dimensional twist-2 operators, the domain of small longitudinal distances is de- 
termined by known non-perturbative QCD methods. However for increasing longitudinal 
distances more and more twist-2 matrix elements are required and soon one enters a region 
where lattice QCD and other non-perturbative approaches are practically not applicable 
any more. Phenomenologically, for unpolarized charge conjugation even quark and gluon 
distributions one observes a transition between the region of small and large longitudinal 
distances around a length scale which corresponds to the electromagnetic nucleon size 



of around 2 fm |13[ . Beyond that scale the matrix elements of the corresponding string 
operators become smooth and approximately flat as functions of longitudinal distance - 
at least at low normalization scales. The matrix elements of the string operators which 
correspond to C-odd unpolarized quark distributions decrease and become small at lon- 
gitudinal distances beyond 2 fm. In both cases, however, the region of large longitudinal 
distances is beyond the scope of presently available non-perturbative methods, and one 
has to resort to approximations such as Regge theory. An ideal non-perturbative QCD 
observable should therefore avoid contributions from large longitudinal distances, i.e. it 
should be sensitive only to those degrees of freedom which reveal themselves at length 
scales which are not larger than the nucleon size. 

Although the total gluon polarization receives most interest from a phenomenological 
point of view, establishing a framework for its evaluation has its own importance. This is 
due to the fact that there is no local, gauge invariant operator which can serve as a "gluon 
polarization partonometer" , i.e. which yields a matrix element associated with AG(/i 2 ). 
In general two different ways to define AG(/i 2 ) have been discussed in the literature so 
far. It has been known for a long time that the gluon polarization can be related to the 
matrix element of a gauge invariant, but non-local gluonic string operator [14, PI : 



NL = n^n u dXTrG^{Xn)[Xn;0]G^{0). (2) 

J 

Here and G^ v denote the gluon field strength and its dual, respectively, and is a 
light-like vector with n 2 = 0, n-a = a° + a 3 = a + for any four-vector a. The trace in (Q) is 
performed in color space and the path-ordered exponential in the adjoint representation 
[An; 0] guarantees gauge invariance. Note however that so far most experience has been 
obtained in dealing with matrix elements of local operators, and therefore it is not clear 
a'priori how one should apply the non-local operator (|2|) in practice. In our recent work 
TH] we have argued that a computation of AG(u 2 ) using the operator Onl at some low 
normalization scale is possible, but requires an insight into the nature of contributions 
arising from large longitudinal distances. Once one accepts a point of view, supported 
e.g. by Regge theory, that despite its non-local character the gluon polarization AG(/i 2 ) 
receives only minor contributions from large longitudinal distances, an approximate "gluon 
polarization partonometer" can be constructed in a gauge-invariant way. In its simplest 
version it takes into account information encoded in only two first computable QCD 
moments of the polarized gluon distribution function. 

On the other hand in light-cone gauge n ■ A = A + = the operator Onl assumes a 
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local form, identical to the n ■ K = K + component of the topological current [II , |T5f : 

= ^ W^a (d^Al + l -gf abc A p b A{) . (3) 

Consequently it is suggestive to use n ■ K as a gluon polarization partonometer due to 
its local character. Corresponding calculations in the framework of the bag model have 
been performed recently in []17|]. Nevertheless, as we shall discuss, although formally the 
operator n ■ K is built from local fields, its matrix element is sensitive to large longitudinal 
distances in the same way as the matrix element of the non-local operator in (||]). Thus 
the advantage of using n ■ K instead of Onl is illusory. This is illustrated below within the 
framework of QCD sum rules. We find that our estimates for AG(fi 2 ) as obtained from 
the gluon polarization partonometer presented in |TB[, and from the local operator (|3|) 
differ approximately by a factor of four. This discrepancy emphasizes, as we will show, 
the role of contributions from different longitudinal length scales. 

The remainder of this paper is organized as follows: in Sec. 2 we collect the most 
important facts about the operator definition of the polarized gluon distribution and the 
gluon polarization integral. In Sec. 3 we present a QCD sum rule estimate for AG(/z 2 ) 
starting out from the matrix element of n ■ K in light-cone gauge, and review an estimate 
for AG(/i 2 ) using the gluon partonometer introduced in Ref . fT6fl . We then discuss and 
explain in Sec. 4 the reasons for the surprising discrepancy between the results obtained 
with these two methods. Finally Sec. 5 is devoted to a summary and conclusions. 



2 Polarized gluon distribution in QCD 

In QCD parton distributions can be related to matrix elements of twist-2 non-local op- 
erators [Q. In this framework unpolarized and polarized gluon distributions are defined 
through matrix elements of the light-cone string operators: 

O G (A;0) = n^n u TrG^(A)[A;0]G^(0), (4) 

O AG (A;0) = rv<TrG^(A)[A;0](V(0). (5) 

Here A stands for a light-like vector being proportional to n. The path-ordered exponen- 
tial 



[A; 0] = P exp 



igAj dXA^(AX) 
Jo 



(6) 



with the strong coupling constant g and the gluon field A 11 guarantees gauge invariance 
of the parton distributions The forward matrix elements of the string operators 

([IH) between nucleon states with momentum p and spin s define the unpolarized and 
polarized gluon distribution of a nucleon, G(-u,/i 2 ) and AG(u, function of the 

Bjorken variable u and the normalization scale /i 2 : 

1 J2(P^\O G (A,0)\p,s)^2 = (p ■ n) 2 du u G(u, fi 2 ) cos \u(p ■ A)] , (7) 
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(Pi s|Oag(A> 0)|p, s)„2 = (p ■ n)(s ■ n) / duu AG(u, fj 2 ) sin \u(p ■ A)] . (8) 

Jo 
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The invariant measure of the light-cone distance between the two gluon fields in is 
given by the so called Ioffe-time z = p ■ A. Furthermore note that for a target polarized 
in the 3-direction one has s ■ n = p ■ n. Taking the Fourier transform of Eqs. (0, ||) yields 
the distribution functions: 



\ roo 

uG(u,/j 2 ) = / dz ^2(p, s\O G (A,0)\p, s) M 2cos( 

7T I p ' f~b J J q 



UZ) 



uAG(u,n 2 ) 



7T (p ■ n)(s ■ n) Jo 



dz (p, s|0 AG (A, 0)|p, s)„2 sia(uz). 



(9) 
(10) 



These definitions are of course in agreement with the perceptions of the parton model. 
Indeed in light-cone gauge, n ■ A — 0, one can express the distributions in terms of 

right- and left-handed gluon operators defined as G +R ^ = e R ^G +fJj , with the polariza- 
(0, -1, -i, 0)/y/2 and el = (0, 1, -i, Q)/y/2: 



tion vectors 



u G(u, /j 2 



dX cos (p ■ n Xu) 



Trp ■ n Jo 

x J2 Tr ^ s \ (G +R (nxW G +R (0) + (G +L (nX)V G+ L (0)\p, s)?, 



'IF 



uAG(u,n 2 



2i 



Ti s ■ n 



roo 

/ dX sin (p ■ n Xu) 
Jo 



x Tr(p,s\ (G + *(nA))V fl (0) - (G +L (nA)) f G +L (0)\p, s)^. 



(12) 



After rewriting the gluon field strength tensor in terms of light-cone quantized fields one 
obtains: 



u G(u, /i 2 ) 
u AG(u,/i 2 ) 



dx d 2 k±5(u — x)x n g (x, k±, R, /x 2 ) + n g (x, k±, L, /i 2 
dx d 2 k±5(u — x)x n g (x, k±, R, /i 2 ) — n g (x, k±, L, fj 2 



(13) 
(14) 



where n g (x, k±, R(L), /i 2 ) denotes the light-cone distribution function of right- (left-) 
handed gluons with light-cone momentum fraction x and transverse momentum kj_ ||14|| . 

In the following we focus on the polarized gluon distribution. Performing a Taylor 
expansion of (|8]) around A = leads to well-known relations, or sum rules, between the 
moments of the polarized gluon distribution AG(u, fi 2 ) and nucleon matrix elements of 
local QCD operators ||18|| : 

f 1 dun 1 ' 1 AG(u,fx 2 ) = TiU 2 ), with I = 3,5,..., 
Jo 

n„n v Tr (p, s\G^(0)(in ■ D) 1 ' 2 ^ »(Q)\p, s) = (s ■ n)(p ■ n) 1 T^ 2 ). (15) 

Note, however, that a sum rule for I = 1 which correspond to the integrated gluon polar- 
ization AG(fi 2 ) is lacking. This is due to the fact that a suitable gauge invariant, charge 
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conjugation even, local operator which may serve as a gluon polarization partonometer 
does not exist. 

On the other hand we find from (|8|) that the polarized gluon distribution AG(fi 2 ) is 
determined by an integral over the corresponding Ioffe-time distribution which is defined 
as: ^ 

T(z,ij 2 )= / du uAG{u, /i 2 ) sin(wz). (16) 
Jo 

Since in the distribution sense one has: 

00 1/1 1 \ 1 

dzsm(uz) = -[ + )=PV-, (17) 

2 \u + ie u — it) u 

where PV denotes the principal value prescription, we indeed obtain: 

AGV) = / dzT(z,fx 2 ). (18) 
J 

It is important to realize that T(z,fi 2 ) ~ z a ~ 2 for large z, if AG(u,/j, 2 ) ~ u~ a at small 
u. Therefore as long as a < 1 the integral over the polarized gluon distribution (|T8| ) 
exists as it converges at large z in an absolute sense. (At small z the integrand in (]l8f) 
should not cause harm since there r(#) ~ zr 3 , with an anticipated finite third moment 
T 3 .) Expanding the RHS of fll6|) around 2 = yields a Taylor expansion of the Ioffe 
distribution T(z,^ 2 ) with coefficients proportional to the odd moments T;(/i 2 ): 

T(z, v 2 ) = T 3 ^ 2 ) z - lr 5 Gu 2 ) + y^Gu 2 ) 2 5 - . . . (19) 

Since each of these moments can be calculated, at least formally, as a reduced matrix ele- 
ment of a local, gauge invariant operator, the convergent integral dl8f) and hence AG(/i 2 ) 
itself is a gauge-invariant quantity. 

As already mentioned in the introduction, in light-cone gauge the gluon polarization 
AG(fi 2 ) can also be related to the expectation value of the topological current (0): 

(p,s\n-K(0)\p,s)^ = ( a .n)-i / dzT(z^ 2 ) = (s ■ n)-^AG(/i 2 ). (20) 

27T JO 27T 

Indeed, in light-cone gauge the operator Oag(A; 0) in (J|) reduces to: 

Oag(A; 0) = ^A 2 (A)9 + A 1 (0) - a+A 1 (A)a+A 2 (0) , (21) 

where A 1 and A 2 denote the transverse components of the gluon field. This simplification 
occurs because in this gauge the path-ordered exponential [A; 0] = 1, and the components 
of the gluon field strength tensor which enter Oag(A; 0) are given by G +± = d + A ± . Using 
the definition (IIBD we obtain: 



1 POO 

AGV) = - ( V( r dz (p, s\d + A (A)d + A l (0) — d + A 1 (A)d + A 2 (0)\p, s)^ , 

2{s ■ n)[p ■ n) Jo 

dX(p,s\—A 2 (Xn)d+A 1 (0)-—A 1 (Xn)d + A 2 (0)\ P ,s) ll 2. (22) 



2(s-n) Jo d\ v ' OX 
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Integrating by parts and assuming that the boundary term vanishes for A — > oo, one 
finally obtains: 

■ n)AG(/2 2 ) = (p, 8^(0)^ A 2 (0) - A 2 (0)d + A 1 (0)\p, a), (23) 

which is easily shown to be equivalent to fl20p . This relation involves an operator which 
is built from local fields. Therefore its matrix element is calculable using known methods 
to deal with local operators. This fact however turns out to be of no real advantage as 
we will point out in Sec. 4, where we describe a QCD sum rule calculation. 

To complete our discussion we derive the scale dependence of the gluon polarization 
AG(/i 2 ) starting out from the operator definition ([18|) . The one-loop evolution equation 
for the Ioffe-time distribution T(z,fi 2 ) reads |L9], pD| : 



where: 



^ ff ' = J Q dv[K x {v)T{vz^ 2 ) + K 2 {v)-p{vz^% (24) 



K x {v) = ^8{v) + 2C A Uvv + ^--l-v\, 

K 2 (v) = -C F (8(v)-2v), (25) 

with Cp — 4/3, = 3, fa = 11 — |iV/, with Afy being the number of active flavors, and 
u = 1 — u. The factor 1/z in ( p4]) arises in accordance with the definition of S(z,/i 2 ), the 
Ioffe-time distribution of polarized quarks: 

(p,s|*(A)n T 5[A;0]*(0)b,s> + (A -> -A) = 4(s • n)£(*, /i 2 ), 

with £(z,/i 2 ) = / du Aq(u, ji 2 ) cos(uz) . (26) 
•/ o 

Here Aq(u, fi 2 ) is the flavor-singlet polarized quark distribution in momentum space. The 
evolution equation for AG(/i 2 ) is obtained by integrating both sides of Eq. fl2"4|) over z. 
We find that gluons enter only via the term in K\ being proportional to /3q. Furthermore 
the quark contribution can be transformed conveniently using the identity: 

— ( dv (S(v) — 2v) cos (uvz) = —2u [ dvv (l v) sin (uvz). (27) 

z Jo Jo \ 2 / 

Applying Eq. (|i"7|) yields then the standard one-loop evolution equation |l| 



47T 

with i = log(/x 2 ), and a s = — — . 

p t 



Here 



S(^ 2 ) = / duAq(u^ 2 ) (29) 



o 
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is the quark polarization in the proton]]. At the one loop level one has fl~8 

dS(// 2 ^ 



dt 



0, and thus ) = S = constant, (30) 



i.e. to this accuracy the quark polarization is scale-independent. One then obtains as a 
solution of <M) the well known result 1T8|, El: 



AGV 2 ) = AGV) + -So - 1 . (31) 

Thus we have shown that the operator definition of AG(/i 2 ) is indeed equivalent to parton 
model considerations. 



3 Different gluon polarization partonometers 

As discussed above, the gluon polarization in the nucleon can be obtained either from 
the integral over the Ioffe-time distribution (|18|), or in light-cone gauge from the for- 
ward matrix element of the topological current (|20|). Here we explore and compare both 
approaches in the framework of QCD sum rules. 



3.1 Gluon polarization from the topological current 

To estimate the nucleon matrix element of the topological current (^) we perform a stan- 
dard QCD sum rule calculation. For this purpose we consider the three-point correlation 
function 

I K = i 2 j d A xe^ x J d^e^iOlTi^n- K(y)fj G (0))\0) (32) 

of the operator n • K in light-cone gauge and the nucleon interpolating currents r/ G , fj G . 
For the latter we take: 

2 

Vg(x) = ^(Void(x) -T/exW), (33) 
Vom(x) = e abc (u aT (x)C lfM u b (x)) l5l ^ al 3 [gG^(x)d(x) 
VeK {x) = e abc (u aT (x)C^d\x))^a aP [gG afi (x)u(x) 

Its overlap with the nucleon state, 

(Q\VG(0)\p,s) = m 2 N \ G u(p,s), (34) 

at the scale /i 2 ~ 1 GeV 2 has been determined in Ref . ||22|| . Note that this current, which 
contains explicit gluon degrees of freedom, has been successfully employed in investigations 

3 For considerations of the one-loop evolution equation for AG(/i 2 ) it is not necessary to discuss the 
role of the axial anomaly in the interpretation of £(/i 2 ). 
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of nucleon matrix elements of QCD operators being sensitive to gluon components of the 
nucleon wave function [^, |23], j24| . To avoid large t-channel contributions we stay at 
Euclidean momenta Q 2 = —q 2 ~ (1—4) GeV 2 and perform a numerical extrapolation 
to Q 2 = at the end. Furthermore the kinematic is chosen such that q ■ n — 0, i.e. 
q 2 = -q ± 2 . 

In the following we concentrate on the contribution to the correlator (|32|) of the form:[] 

l,p{p-n)T K {p 2 ,{p + q) 2 ,Q 2 ). (35) 

The invariant function T% can be projected out uniquely from It receives contributions 
from nucleons as well as from higher resonances and continuum states. For the nucleon 
contribution we have: 

T K (P ,(p + q) ,Q) = = {m% _ f)im% _ {p + q)2) ^G(Q ), (36) 

where the form factor AG(Q 2 ) coincides at Q 2 = with the gluon polarization AG. 



As a next step we use the fact that admits a double spectral representation [Z5(j: 

T K (P 2 , (p + q)\ Q 2 ) = I ~^~ 2 I p— 2 Pk(su 82, Q 2 ). (37) 

J s 1 -p 2 J s 2 — [p + q) 2 

We have calculated the spectral density p K in light-cone gauge, taking into account the 
dimension-1 and dimension-6 operators in the OPE of Tk, which are visualized in Fig. 1. 
We obtain: 



m o Oi s Q 2 (Aa - 

Pk\s u s 2 ,Q*) = * VV 6 ' (A q Q 2 Rl-3A q R 2 -Q 2 R 2 

23040 7T b Rz 



.3 



- 7R* -2A q Q z Rz s l s 2 + 2A q Rs 1 s 2 + AR*s l s 2 -AQ l s{s z 2 
pV(si)«2,Q) = -5 — W) -s , (38) 

9 7T Rl 

where A q = si + s 2 + Q 2 and i? = A 2 — 4sis 2 . To eliminate contributions from higher 
resonances and the continuum we limit the integral over s\ and s 2 by the continuum 
threshold sq. Performing in addition a Borel transformation in p 2 and (p + q) 2 yields the 
sum rule: 

^AG(Q 2 ) = e - T -^- / ds x / ds 2 e-^y™ 2 p K (s h s 2 , Q 2 ) , (39) 

Z7T *G m N ~' ~' 

where M 2 is the Borel parameter. For consistency both, s and M 2 should be taken around 
their values fixed by the two-point sum rules p2"]. In the actual calculation we use at the 



4 We suppress here and in the following any dependence on the scale /i 2 
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scale /i 2 = 1 GeV 2 the standard value for the quark condensate, — (27r) 2 (gg) = 0.67 GeV 3 , 
and the strong coupling a s = 0.37 f22fl . 

The stability of the sum rule ( |39"D against variations of the Borel parameter M 2 and 
the continuum threshold sq is illustrated in Figs. 2 and 3. One observes a relatively strong 
dependence on sq which can be traced back to the large dimension of the interpolating 
current rjc- To determine the gluon polarization AG we extrapolate the form factor 
AG(Q 2 ) to Q 2 -»■ 0. For this purpose we fit the RHS of in the interval 1 GeV 2 < 
Q 2 < 4 GeV 2 by: 

a<5(<?2) = ag ^TqW (40) 

with 7 = 3 as suggested by quark counting rules. Extrapolating to Q 2 — > gives AG{ji 2 ~ 
1 GeV 2 ) = 0.6 ± 0.2, where the error has been estimated from the M 2 and sq dependence 
of the extrapolation. We have checked that this result essentially does not change if we 
allow 7 as a fit parameter as well. An additional 30% error arises from the uncertainty in 



Xc and the vacuum saturation ansatz for the four-quark condensate [22 



3.2 Gluon polarization from an approximate Ioffe-time distri- 
bution 

In comparison to the above calculation we review here an estimate of AG via an ap- 
proximate Ioffe-time distribution [|l(J. The latter is based on the conjecture that in the 
laboratory frame the polarized gluon distribution (|I3|) receives major contributions only 
from longitudinal distances smaller than the nucleon diameter, as determined by the 
nucleon electromagnetic form factor. This hypothesis is supported both by Regge phe- 
nomenology and the color coherence hypothesis [p7 |, which impose strong restrictions 



on polarized glue at small u or, equivalently, at large z. For a more accurate determination 
of the longitudinal length scale at which the gluon polarization ( |T8| ) effectively saturates 
an additional assumption has to be made. In this respect we have constructed in a 
polarized gluon correlation function T(z) which behaves at large z similar to the valence 
quark distribution of a nucleon, i.e. it becomes small for z > 10. (In the nucleon rest 
frame z — 10 corresponds to a longitudinal distance of 2fm.) Indeed most of the current 



parametrizations for polarized glue sustain this picture f28fl , as shown in Fig. 4 for the 
distributions from Refs. [P7|, E^|. 

In Fig. 5 we present the Ioffe-time distribution as obtained from the expansion in 
Eq. (|i9|) , taking into account terms up to the Z-th moment of AG(u). One finds that the 
first two non- vanishing moments T3 and determine the Ioffe-time distribution T(z) at 
small z nearly up to its maximum^. Therefore if the polarized gluon distribution is of 
regular shape similar to the ones shown in Fig. 4, and obtains significant contributions 
only from the region z < 10, a simple estimate of AG(/i 2 ) is possible [|16 . It is given by 



the area of the triangle spanned by the points z = 0, z = 10, and the maximum of the 



5 Note that in the notation of Um T 3 and T 5 correspond to T 2 and T 4 , respectively. 
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approximate Ioffe-time distribution 



T(z) « T»z 



6 



r, z\ 



(41) 



This rough estimate requires only the knowledge of two moments of AG(u). These can, 
as a matter of principle, be taken from any theoretical investigation, i.e. QCD sum rules 
or lattice calculations. 

Of course one can do better if the normalization of T(z) is known at some large value 
of z ~ 10. Then a more accurate estimate can be achieved, calculating the area bound 
by the approximate Ioffe-time distribution ([4i~D up to its maximum and a straight line 
connecting this point with the value of T(z = 10). For the model parametrizations [27, 



this leads to an estimate for AG with (10 — 20)% accuracy. In [16| the color coherence 
hypothesis |27J has been applied to obtain information on the large z or, equivalently, 
small u behavior of polarized glue in nucleons: 



AG(u) 
G{u) 



u. 



for u — > 0. 



(42) 



It allows to estimate T(z = 10) from parametrizations for the unpolarized gluon dis- 
tribution. For example one finds from the GRV [ 30 1 and CTEQ |31| LO unpolarized 
gluon distributions T(z = 10, /i 2 ~ 1 GeV 2 ) = 0.005 - 0.007. On the other hand the 
parametrizations of |7j and |29[ yield T(z = 10, /j 2 ~ 1 GeV 2 ) = 0.02 and 0.09, re- 
spectively. Since in principle nothing prevents T(z) from becoming negative at large z, 
—0.05 < T(z = 10) < 0.05 should be a conservative estimate. 

The moments T3 and have been calculated within a standard QCD sum rule ap- 
proach, starting from an investigation of the three-point correlation function (TBI: 



Ir = t 2 Jd A xe^ x Jd'ye^yiOlTi^OAGiy + A^y-A^GmiO), 

= 1 ,p(p-n) 2 T T (p 2 ,(p + q) 2 ,Q 2 ,z) + .... (43) 



Separating the contribution from nucleon states yields: 

Tr {nj 5 I T ) 



T T (p\(p + q)\Q\z 



X 2 G m% 



4(p • n) 3 {m 2 N - p 2 )(m 2 N - (p + g) 2 ) 



f(^,Q 2 ) + 



(44) 



where the form factor T(z,Q 2 ) coincides at Q 2 = with the Ioffe-time distribution T(z) 
in ([T5|) . In a next step again a double spectral representation is introduced: 



Tr(p 2 , (P + d)\ Q\ z)=f -^- 2 j p— 2 P r(su Q\ z). 

J Si - p l J s 2 — (p + q) 1 



(45) 



As in Sec. 3.1 the dimension-1 and dimension-6 contributions have been considered. For 
them the polarized spectral density pr coincides with the spectral density for the operator 
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(||) which determines the unpolarized gluon distribution G(w,/i 2 ) p2"| . Since this spec- 
tral density leads to a realistic large value for the contribution of gluons to the nucleon 
momentum, it is suggestive that also the integrated gluon polarization is large. 

14|,[45D in powers of z yields sum rules for ^(Q 2 ) and T^(Q 2 ), which 
]. An extrapolation to Q 2 = as in ( |4"0"D finally leads to the de- 
sired moments. In combination with the color coherence hypothesis (|42"D one then obtains 
AG(/i 2 ~ 1 GeV 2 ) = 2 ± 1. The main source of the quoted error is due to uncertainties 
in the QCD sum rule approach. Furthermore, under the assumption that T(z) follows 
the Regge behavior at large z, possible contributions from longitudinal distances beyond 
the ones accounted for have been estimated to be smaller than ±0.2. Note however that 



Expanding Eqs. 
are given explicitly in 



current experience p2j, |33[ indicates that QCD sum rules significantly overestimate higher 
moments of parton distribution functions. Therefore the result quoted above should be 
treated as an upper limit for AG(/i 2 ~ 1 GeV 2 ). 



4 Discussion 

We have found that the gluon polarization obtained from the matrix element of the 
topological current differs by a factor of around four as compared to the estimate based 
on the approximate Ioffe-time distribution. This obviously calls for an explanation. As we 
show below, this discrepancy emphasizes the role played by contributions from different 
longitudinal distances. 

First let us comment on the question whether we have gained anything by going from 
the non-local definition of AG to the local one. To this end note that due to Eq. (|20|) one 

Pk(si,s 2} Q 2 ) = dz p r (s 1 ,s 2} Q 2 } z) . (46) 
Jo 

Of course this relation can be directly verified for the contributions from the operators of 
dimension one and six as used in Sec. 3.1 and 3.2. This however implies that, in accordance 
with the equivalence of the definitions ([18]) and (p0|), the matrix element of n ■ K is as 
sensitive to different longitudinal distances as the matrix element of the non-local operator 
([5]). Consequently the apparent local form of K does not remove the essential non-local 
character of AG. 

To clarify the reason for our different results for AG consider a polarized gluon distri- 
bution with the following behavior: 

( Au~ a , for it — > 0, 

G(u) ~ I (47) 
[ B{l-uf, for u-> 1. 

The asymptotic form of the corresponding Ioffe-time distribution at large z arising from 

6 We remind the reader that this relation holds only if the gauge-dependent density px is calculated 
in light-cone gauge. 
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the regions u — > and u — > 1 reads: 

r(z) = -4sin (^a)r E (2 - a)z a ~ 2 + ...-Bcoa(z- /30r E (/3 + l)^ 1 ^ + . . . , (48) 

where denotes the Euler Gamma function. Phenomenologically one would expect 
< a < 1 and (3 > 4 ||27|| , which guarantees a smooth non-oscillatory behavior of T(z) at 
large distances. (In fl48|) even for a ~ the next-to-leading contribution from the small 
u region dominates over the leading one arising from the region u — > 1, if only (3 is large 
enough and A and B are of similar magnitude.) 

The Ioffe-time distribution obtained from QCD sum rules behaves however differently. 
This can be seen from Fig. 6 where we show its dominant contribution which results 
from the dimension-6 operator taken in the limit Q 2 — > 0. This limit should give a 
qualitatively reasonable estimate since the Q 2 -dependence of the sum rules in Sec. [Tl 
and |3.2j has turned out to be quite smooth. We find that, contrary to what one expects 
from phenomenological considerations, T(z) oscillates strongly and decreases relatively 



slow at large z. In terms of the characteristic exponents in (|48|) its behavior corresponds 
approximately to a ~ — 1 and (3 ~ 0. This shows that QCD sum rules yield a distribution 
skewed towards large values of u, and consequently leads to predictions for the moments of 
AG(u) which are too large. The small value of (3 reflects a large weight for configurations in 
which one gluon carries most of the momentum of the nucleon. In the present calculation 
this is due to the fact that neither a perturbative (Sudakov), nor a non-perturbative (large 
invariant mass) suppression mechanism for such configurations is present in lowest order 
OPE. Since the sum rule for n ■ K in Sec. p.l| accounts for the full integral over T(z) from 
zero to infinity, it incorporates all oscillations at intermediate and large z, leading to a 
small value of AG. Apart from the small- z region the z-dependence of T(z) shown in 
Fig. 6 is certainly not realistic. Therefore we have to conclude that the QCD sum rule 
calculation of the matrix element of n ■ K is not entirely self-consistent, as it receives large 
contributions from regions where the sum rule method is not reliable. 

This is different in the approach discussed in Sec. |3.2| . Here we have assumed that 
QCD sum rules yield a reasonable estimate for the first two moments T3 and Y§. However 
we have discarded oscillations of T(z) at large z, assuming that the main contribution to 
AG arises from small light-cone distances. As a consequence QCD sum rules are used 
only in a domain where they are in principle applicable. We therefore believe that such 
an estimate of AG is better justified. 



5 Summary and conclusions 

In QCD the twist-2 polarized gluon distribution is defined through a gauge invariant 
but non-local string operator. As a consequence AG can receive, at least in principle, 
contributions from different longitudinal distances. Although in light-cone gauge AG can 
be formally expressed through the forward matrix element of the local topological current, 
also in this case contributions from all longitudinal distances are accumulated. Therefore 
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the latter can be used for a trustworthy estimate of AG only if an approximation to 
strong interaction dynamics is available which is applicable at all longitudinal distances. 
We have illustrated this point in the framework of a QCD sum rule calculation which 
leads to rather unrealistic contributions from large longitudinal distances, and results in 
a small value for AG(/i 2 ~ 1 GeV 2 ) = 0.6 ± 0.2. 

For a self consistent estimate of the gluon polarization in nucleons one has to ensure 
that the main support to AG results from distances where the used approximation is sup- 
posed to do its best. In the case of QCD sum rules only contributions from small distances 
can be approximated in a reasonable way. Combining the latter with the assumption that 
contributions to AG from distances larger than the typical nucleon size are small yields 
AG(/i 2 ~ lGeV 2 ) = 2± 1. 

If, contrary to our assumption, AG receives at low normalization scales fi 2 important 
contributions from large longitudinal distances, new interesting questions would arise since 
such contributions could hardly be interpreted as being due to confining glue, understood 
as part of a nucleon with a size of around one fm. At large normalization scales GLAP 
evolution may result in a more and more singular behavior of AG(u,/i 2 ) at small values 
of u. In this case it is understood, however, that such contributions have to be filtered 
out in order to learn something about the distribution of nucleon polarization among 
low-virtuality degrees of freedom. 

To avoid in the calculation of moments of a parton distribution a situation in which 
an approximation tailored for contributions from small longitudinal distances generates 
large contributions from distances beyond its scope of applicability, one should require 
that it gives a reasonable behavior of the considered distribution at large values of z or, 
equivalently, at values of u close to 0[] If such a requirement cannot be fulfilled an explicit 
construction of contributions from large longitudinal distances , like in []16||, is necessary. 
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7 To achieve this for example for the quark distribution in the bag model one has to introduce a com- 
plicated projection which allows to approximately express bag model solutions in terms of eigenfunctions 
of the momentum operator pm . 
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Fig. 1: Graphs representing the contributions from dimension- 1 and dimension-6 opera- 
tors to the QCD sum rule calculations described in the text. 
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2: Stability of the sum rule for the nucleon matrix element of n-K against variations 
of the Borel mass between M 2 = 1 GeV 2 (lower curve) and M 2 = 2 GeV 2 (upper 
curve). The continuum threshold has been fixed at y/s = 1.5 GeV. The dots 
correspond to the sum rule results (|39D for 1 GeV 2 < Q 2 < 4 GeV 2 . The solid lines 
show the extrapolation to Q 2 — > (^0|). 
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Fig. 3: Stability of the sum rule for the nucleon matrix element of n ■ K against varia- 
tions of the continuum threshold between yfs Q = 1.4 GeV (lower curve) and \fs = 
1.6 GeV (upper curve). The Borel mass has been fixed at M 2 = 1.5 GeV 2 . The dots 
correspond to the sum rule results (|39"D for 1 GeV 2 < Q 2 < 4 GeV 2 . The solid lines 
show the extrapolation to Q 2 — > (flop. 
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Fig. 5: The Ioffe-time distribution of polarized glue, approximated by its n-th order 
Taylor expansion around z = (|I9|) , using the parameterization of Ref.[|29 . 
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